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ABSTRACT 



In this first paper on the results of our Herschel PACS survey of local Ultraluminous Infrared Galaxies (ULIRGs), as part of our 
SHINING survey of local galaxies, we present far-infrared spectroscopy of Mrk 231, the most luminous of the local ULIRGs, and a 
type 1 broad absorption line AGN. For the first time in a ULIRG, all observed far-infrared fine-structure lines in the PACS range were 
detected and all were found to be deficient relative to the far infrared luminosity by 1 - 2 orders of magnitude compared with lower 
luminosity galaxies. The deficits are similar to those for the mid-infrared lines, with the most deficient lines showing high ionization 
potentials. Aged starbursts may account for part of the deficits, but partial covering of the highest excitation AGN powered regions 
may explain the remaining line deficits. A massive molecular outflow, discovered in OH and '*0H, showing outflow velocities out to 
at least 1400 km sec"' , is a unique signature of the clearing out of the molecular disk that formed by dissipative collapse during the 
merger. The outflow is characterized by extremely high ratios of '**0 / "'O suggestive of interstellar medium processing by advanced 
starbursts. 
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1. Introduction 

ULIRGs (L/s > IO'^Lq), which in the local universe are known 
to be mergers of gas rich galaxies, are thought to play a critical 
role in galaxy evolution. It is therefore crucial to understand the 
conditions, dynamics, chemistry and energetics of this stage of 
evolution, in which dissipative collapse accompanies the trans- 
formation of gas rich galaxies into ellipticals (Lonsdale et al. 
2006). As part of the SHINING Key Project on local galaxies, 
we are carrying out a far-infrared (FIR) spectroscopic survey of 
all 21 ULIRGs in the Revised Bright Galaxy Survey (Sanders 
et al. 2003). Previous work on ULIRGs with the Infrared Space 
Observatory (ISO) Long Wavelength Spectrometer (LWS) found 
deficits in far-infrared atomic and ionized fine-structure line 
emission relative to their FIR luminosities (Luhman et al. 1 19981 1 
accompanied by prominent molecular absorption in excited tran- 
sitions of molecules such as OH and H2O rarely seen in other 
galaxies (Fischer et al. I1999I I. The deficits, often based only 
on a few non-detections, were postulated to result from the 
high ratios of UV radiation density to particle density in the 
nuclei of these galaxies (MaUiotra et al. 2001, Luhman et al. 
2003, Abel et al. 2009), high gas density (e.g. Negishi et al. 
2001) and high FIR opacity and/or high luminosity-to-mass ratio 



* Herschel is an ESA space observatory with science instruments 
provided by Principal Investigator consortia. It is open for proposals 
for observing time from the worldwide astronomical community. 
** Visiting scientist at Max-Planck-Institut fiir extraterrestrische 
Physik (MPE), Garching, Germany 



(Gonzalez-Alfonso et al. 2004, 2008, hereafter GA04, GA08), 
but because of the paucity of diagnostic line detections this puz- 
zle remained unresolved. GA04 and GA08 showed that the high 
excitation molecular absorptions in both Arp 220 and Mrk 231 
are radiatively pumped, consistent with high radiation density to 
particle density. They identified far-infrared absorption by OH, 
H2O, NH3, NH, and '^OH, but the identifications of the latter 
species were uncertain due to the low resolution of the LWS. 
The PACS ULIRG spectroscopic survey is designed to kinemat- 
ically identify the ionized, atomic and molecular regions and to 
study the conditions of the interstellar medium, thereby illumi- 
nating the nature of this important evolutionary phase. In this 
first paper we present our initial results on Mrk 231, the most 
luminous of the local ULIRGs (L(8-1000A(m) = 3.2 x lO'^ Lq) 
and a type 1, low-ionization broad absorption line (LoBAL) ac- 
tive galactic nucleus (AGN) at an adopted distance of 172 Mpc 
(z=0.04217). Its central quasar is covered by a semi-transparent 
dusty shroud producing about 3.1 magnitudes of extinction at 
4400 A (Reynolds et al. 2009) and is at the center of a rotat- 
ing, nearly face-on molecular disk (Downes & Solomon 1998). 
Based primarily on Spitzer results, Veilleux et al. (2009) estimate 
that the average AGN contribution to the bolometiic luminosity 
in ULIRGs is 35 - 40% and that for Mrk 23 1 the AGN contribu- 
tion is ~ 70% by most estimation techniques. The contribution 
of an advanced 120 - 250 Myr nuclear starburst is estimated at 
25 - 40% based on near infrared observations of Mrk 23 1 by 
Davies et al. (2007). 
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Table 1. Spectroscopic measurements for Mrk 23 1 



Transition"- * 


Line flux 


Line flux en'or^ 


Equivalent Width 


Measured FWHM 


Inferred FWHM" 




10-" Wm-2 


10-"Wm-2 


10-3 pm 


kms-' 


kms-' 


[N III] 57.3 //m 


2.8 


0.6 


-1.0 


203 


177 


[01] 63.2 //m 


36.0 


2.6 


-17.0 


232 


218 


[0 III] 88.4 //m 


4.1 


0.6 


-3.8 


311 


289 


[Nil] 121.9 pm 


4.1 


0.5 


-10.4 


389 


266 


[O I] 145.5 nm 


3.2 


0.4 


-17.9 


324 


208 


[CII] 157.7 pm 


38.3 


1.3 


-296.0 


337 


247 


[Nil] 205.2pm'' 


2.6 


0.2 


-61.0 


524 


242 


OH 79.1, 79.2 pm 


-22.8 


1.0 


15.3 


- 


- 




3.2 


0.6 


-2.1 


— 


- 


OH 119.2, 119.4 pm 


-47.6 


0.3 


107.4 


- 


- 




9.4 


0.3 


-21.6 


- 


- 


'**0H 120.0, 120.2 pm 


-14.7 


0.2 


34.3 


- 


- 




1.7 


0.2 


-4.0 


- 


- 


H2O 78.7 pm 


-18.0 


0.1 


12.2 


444 


422 


HF/H2O 121.7 pm 


-4.7 


0.4 


11.8 


560 


483 



'"* For doublets, fluxes listed are for the blended pair. *** For transitions with both absorption and emission, the absorption component is listed 

first. *''' Errors listed are the statistical rms uncertainties; we allow 25% calibration uncertainties (see text). *''' Inferred velocity width is based on 

the instrument resolution assuming a Gaussian profile; FWHMs are not listed for doublets. *"' Herschel SPIRE data (VIO). 



2. Observations and data reduction 

The observations were taken with the Photodetector Array 
Camera and Spectrometer (PACS) integral field spectrometer 
(Poglitsch et al. 2010) on board the Herschel Space Observatory 
(Pilbratt et al. 2010) in high spectral sampling range spec- 
troscopy mode using small chop-nod cycles. Scans covering the 
range ± 1300 km s-' were taken around five fine-structure lines 
and a longer range scan included the OH 119 pm ^113/2 5/2 - 3/2 
A-doublet transitions, the '**OH 120 pm counterparts, and the 
[Nil] 122 pm fine-structure line. The OH 79 pm ^Hj/i - ^03/2 
1/2 - 3/2 A-doublet transitions, the nearly superposed '^OH 
counterparts, and the H2O 78.7 pm 423 - 3 12 line were observed 
in parallel on the blue aiTay. 

The basic data reduction was done using the standard PACS 
reduction and calibration pipeline (ipipe) included in HIPE 2.0 
134CQ The data are consistent with an unresolved point source. 
For a well centered point source, the peak line and continuum 
flux will fall on the central of the PACS 5x5 spaxels (9" x 9" 
spatial pixels), with nearly all of the rest of the flux falling on 
the remaining pixels of the aiTay. The best signal-to-noise on 
a faint line is obtained by scaling the line profile on the central 
spaxel, to the continuum flux integrated over the array. Where 
better baselines were obtained by scaling the sum of the brightest 
spaxels, this was done instead. For the [CII] 158 pm scan, all of 
the spaxels were co-added. 

At the time of preparation of this paper the full planned cal- 
ibration procedure is not yet incorporated into the pipeline. For 
calibration of the PACS line fluxes, the line scan continua were 
therefore scaled to a fit to existing photometry of Mrk 231. The 
integrated line fluxes, statistical eiTors, equivalent widths, mea- 
sured Gaussian-fit Av/7iy//M and the inferred intrinsic Avp^uM 
obtained by subtracting the instrumental /S.\fWHM in quadrature 
are listed in Table[T] We note that the [CII], [OI]63, and OHl 19 
line fluxes are in agreement with ISO LWS line fluxes reported 
in Luhman et al. (2003) and GA08 to within the statistical un- 
certainties, but the stronger [Nil] line flux reported in GA08 is 
not consistent with our observations. Additionally the [N II] 205 




' HIPE is a joint development by the Herschel Science Ground 
Segment Consortium, consisting of ESA, the NASA Herschel Science 
Center, and the HIFI, PACS and SPIRE consortia. 



-1000 1000 

Velocity (km s ) 

Fig. 1. Mrk 231 normalized far-infrared fine-structure line pro- 
files versus velocity. The red bars represent the instrumen- 
tal FWHMs. The absorption observed to the blue side of the 
[Nil] 122 pm line is due to the HF 2-1 and/or H2O 432 - 423 
lines. 

pm line measured by Herschel SPIRE is included in Tableflland 
Fig. [T] based on measurements presented in van der Weif et al. 
(2010, hereafter VIO). We estimate the calibration uncertainties 

to be ±25%. 



3. Results and Discussion 

3.1. The fine-structure lines 

The continuum normalized fine-structure line spectra are dis- 
played in Fig.[T]as a function of velocity, where Mrk 23rs red- 
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Fig. 2. Line-to-FIR ratio in Mrk 23 1 divided by the median line- 
to-FIR ratio in a sample of HII galaxies (black symbols) and 
AGN (red symbols) versus wavelength {top) and ionization po- 
tential {bottom). The mid-infrared line fluxes in Mrk 231 are 
from Armus et al. (2007). The error bars include the Mrk 231 
calibration uncertainties and the statistical errors for the com- 
parison sample values. 



shift of 0.04217 is defined as km s"' and the red horizontal 
bars represent the instrumental FWHMs. The FWHM velocities 
observed in these lines are in all cases larger than the instrumen- 
tal resolution and we do not detect self-absorption. The average 
of the derived intrinsic FWHMs for the fine-structure lines is 
235 km s"', similar to the FWHM velocity widths of 167 km s"^ 
and 270 km s"' measured by Sanders et al. (1991) and Tacconi 
et al. (2002) from CO(l-O) and near-infrared stellar CO band- 
head studies, respectively. Thus the fine-structure line emission 
is probably associated with the central star forming disk. In the 
[CII] line, wings not characteristic of the instrumental profile are 
present, out to -1000 km s"' for the blue wing, suggesting that 
some C^ is associated with the outflow discussed in the next sub- 
section. The [NII]205 profile also suggests a wing out to -700 
km s"', consistent with the detection of optical [Nil] lines by 
Rupke et al. (2005). We estimate that the absorption by HF or 
H2O observed in the blue wing of the [Nil] 122 //m line does not 
aff'ect the [Nil] line flux by more than 50%, since the FWHM 
of this transition is similar to those of the other fine-structure 
lines. The observed [Nn]205 to [Nil] 122 ratio of 0.6 ± 0.4 trans- 
lates into a density range for N^ of log n^, = 1 .4^g4 (Rubin et al. 
1993). Photodissociation region (PDR) modehng of [CII], [OI], 
and FIR fluxes yields UV radiation field and density Go ~ 6 x 
10^ and n ^ 500 cm"-', at the extreme range of average radiation 
density per particle density seen in extragalactic nuclei (Sturm 
et al. 2010, Fig. 4). 

The fine-structure line strengths relative to the far-infrared 
luminosities are weak compared with normal HII galaxies and 
AGN. In Fig. |2] we plot as a function of wavelength (top) 
and ionization potential (bottom) the fine-structure line to far- 
infrared flux ratio (FIR, 42 - 122 //m) in Mrk 231 relative to the 
median value for subsamples of HII galaxies (black symbols) 
and AGN (red symbols) with luminosities 1 x 10^ Lq < L/« < 
5 X 10" Lq and major isophotal diameters D25 < 250 arcsecs 
compiled from the literature and our SHINING sample (Gracia- 
Carpio et al., in prep.). There are no [Nil] 205 fiva line measure- 
ments for the comparison samples, so this line is not included 
in Fig.|2] The observed deficits are severe: typically an order of 
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Fig. 3. P-Cygni profiles (solid blue, black histograms) of the 
OH 119 yum and '**OH 120 /im partially resolved doublets {top) 
and the OH 79 jum unresolved doublet {bottom) in Mrk 23 1 are 
compared with our modeled profiles (solid red curves, see text). 
Absorption in the H2O 423 - 3 12 transition with a possible blue- 
shifted wing is modeled in GAIO. The zero velocity rest-frame 
wavelengths are marked with black (OH), blue ('^OH), and ma- 
genta (H2O) arrows and the instrumental FWHM is indicated in 
each panel with horizontal red bars. 



magnitude for the neutral and low ionization species and up to 
two orders of magnitude for the higher ionization lines. There 
is no correlation between the deficit and wavelength, and only 
weak correlation with critical density, so differential extinction 
(with wavelength) does not appear to play a significant role, and 
density does not play a major role. We find a strong inverse cor- 
relation (correlation coefficient, -0.78) with ionization potential 
when compared with the AGN sample. A similar inverse corre- 
lation is found with the starburst comparison sample only if the 
[NeIII]15 deficit is treated as an outlier. Low excitation is con- 
sistent with the estimation by Davies et al. (2007) that a 120 - 
250 Myr starburst supplies 25 - 40 % of the bolometric lumi- 
nosity of Mrk 231. However such an aging starburst can create 
a deficit of 40% at most. The high relative strength of the [OI] 
145 yum line, with its high excitation lower level, is suggestive of 
increased radiation density per particle density, as parameterized 
by the ionization parameter U (Luhman et al. 2003; Abel et al. 
2009). The high ratio [OI] 145 / [CII] 158 ~ 1/1 1 lies between 
the predicted line ratio for starburst and AGN at the far-infrared 
F(60)/F(100) color of Mrk 231 with U ^ lO'- (high U) in the 
dust-bounded models of Abel et al. (2009). The strong inverse 
correlation relative to the AGN comparison sample in Fig. [2] is 
consistent with partial covering, higher for high ionization lines, 
of the line emitting region, due for example to a face-on molec- 
ular torus, with optically thick clumps. 

3.2. The molecular outflow 

Figure p] displays the observed 119 - 120 yum OH / '^OH and 
79 yum H2O / OH fine profiles in Mrk 231. The OH lines show 
spectacular P-Cygni profiles in all three ground-state doublets, 
with extremely broad blue-shifted absorption as far out as 1400 
km s"' in the OH 119 fim line. Both the high molecular out- 
flow velocities and the relative strengths of the OH and '^OH 
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are unprecedented. The excited H2O 78.7 //m and HF/H2O 121.7 
yum lines (see Fig. [Til appear to have blue-shifted wings, indi- 
cating that H2O and HP may also be involved in the outflow. 
Recently a molecular outflow was detected by Sakamoto et al. 
(2009) in submillimeter lines of HCO^ and CO toward the nu- 
clei of Arp220 with an outflow velocity of ~ 100 km s"', but this 
is modest compared with that of Mrk 231. 

Our model of the observed profiles is shown in Fig. [3] and 
is based on the continuum components and radiative pumping 
described in Gonzalez-Alfonso et al. (2010, hereafter GAIO). It 
uses three continuum components, hot (radius - 23 pc, T^/,,,, - 
400 - 150 K), warm (120 pc, 95 K), and extended (610 pc, 41 K), 
and three velocity components at 100, 400, and 600 km s"' (the 
latter component with high turbulent velocity, 300 km s"'). We 
have conservatively assumed a screen approach for OH, with the 
outflow located just around the warm continuum source, result- 
ing inA^(OH) - 1.4x10'^ cm"^. The model indicates mechanical 
energy > 10^^ erg (dominated by the high velocity component), 
mechanical luminosity > 1 % of the TIR luminosity, and an out- 
flow mass of > 7 X 10^ M©. Observations of high-lying OH lines 
would help better constrain the outflow location, energetics, and 
the apparent extremely low '^O / "^O ratio. 

Figure0(top) shows that the absorption at 120 ;um by '^OH 
at low outflow velocities (> -250 km s"') is deeper than the 
corresponding absorption in the cross-ladder 79 yum OH dou- 
blet at the same velocity. In this simple model we assume 
that the incomplete absorption in the optically thick OH 119 
//m doublet is a result of the fact that OH at a given veloc- 
ity is only covering a fraction of the continuum (Fig. HI mid- 
dle) and we ignore any differential re-emission of the "*0H and 
79 jum OH lines. We infer that at velocities > -250 km s ' 
T,,(''*OH 120//m) > T,,(OH79//m) (Fig. [ij bottom). From the ra- 
tio of the Einstein B coefficients, we find that the opacity in the 
79 jt/m OH line is 40 times lower than in the 119 /im OH line 
and thus that A^(OH)/A^('^OH) < 40. These results ai-e not criti- 
cally dependent on the uncertainties in the covering factor due to 
re-emission (or thermal emission) in the 1 19 //m OH line, which 
would affect both the opacities of the '^OH and 79 yum OH lines 
but not their relationship. We note that the effects of filling in by 
cooler dust would weaken the 120 /im absorption more than the 
79 //m absorption. The rarer isotopologue '^OH would appear 
between the two lines at 1 19 and 120 yum, but it is not possible to 
put a limit on its strength due to the broad width of the detected 
features. 

We note that CH+, HCl, H^^Cl, H™Br, and H^'Br have tran- 
sitions at « 120 //m, but HCl and HBr and their isotopologues 
are not expected to significantly contaminate the observed ab- 
sorption due to the high lower level energies of their 120 yum 
transitions (180-240 K). On the other hand, CH+(3-2) at 1 19.86 
//m could significantly contaminate the blue-shifted component 
of the "^OH doublet, given that CH+(l-0) is detected in Mrk 231 
(VIO). We estimate that contamination by CH^ could decrease 
the implied '^OH abundance by a factor of 1.5 - 2 and the '*0 
/ '^O ratio could be increased to ~ 30 - 40, still significantly 
below the solar system isotopic ratio of 500, the galactic center 
ratio of 250, the starburst galaxy ratio of ~ 150, and very close 
to the enhancement to ~ 50 that is theoretically predicted with 
a top heavy initial mass function in an advanced starburst stage 
(Henkel & Mauersberger 1993). If our simple model is correct, 
the high '**0 / '^O ratio seen here provides stringent constraints 
for stellar nucleosynthesis models. 

Rupke et al. (2005) find that although both AGN and star- 
burst galaxies have massive outflows traced by broad interstel- 
lar Na I D absorption lines, the high velocity (Av> 2000 km 
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Fig. 4. Normalized velocity profiles for the blue-shifted OH and 
"^OH absorption (top). The velocity scale is given relative to the 
bluest component of each OH doublet. The covering factor {mid- 
dle), assumed to be the same for OH and '^OH and derived from 
the profile of the optically thick OH 119 /im line, is used to cal- 
culate the opacities in the 79 /im and 120 //m lines {bottom, see 
text). 



s ') outflows seen almost exclusively in type 1 ULIRGs includ- 
ing Mrk 231, are probably powered by the AGN. This was also 
shown in models by Narayanan et al. (2008). If it is also true 
for the massive and energetic moderate velocity outflow seen 
here for the first time in OH and '*OH, then an AGN driven out- 
flow is in the process of dispersing into the intergalactic medium, 
highly processed material produced by an advanced starburst. 
The molecular outflow shows velocities similar to the kilopar- 
sec scale outflow component of Mrk 231 seen by Rupke et al. 
and appears to be a unique signature of the clearing out of the 
molecular disk that formed with the dissipative collapse during 
the merger 
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